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Abstract. Much evidence has been collected to date which shows that repulsive electron- 
electron interaction can lead to the formation of particle pairs in a one-dimensional random 
energy landscape. The localization length A2 of these pair states is finite, but larger than 
the localization length Ai of the individual particles. After a short review of previous work, 
we present numerical evidence for this effect based on an analysis of the interaction matrix 
elements and an application of the decimation method. The results based on the decima- 
tion method for a two-dimensional disordered medium support the localization-delocalization 
transition of pair states predicted recently. 
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1 Introduction 

Investigations into the interplay of disorder and many-body interactions continue to 
receive a lot of attention mainly due to the persistent current problem^^"^ and the ex- 
perimental discovery of the two-dimensional (2D) metal- insulator-transition (MIT).''"^ 
In order to theoretically study this interplay, in principle with increasing system size 
one has to solve a problem with an exponentially growing number of states in the 
Hilbert space. At present, this can be achieved only for a few particles in ID^'^ and 
very few particles in 2D.^"^^ Shepelyansky^^'^* — following earlier work of Dhorokov^^ 
— proposed another approach looking at the properties of two interacting particles 
(TIP) in a random environment. The TIP Hamiltonian is 

H = — t (|n, m) (n + 1 , m\ + \n, m) {n, m + 1\ + h.c.) 



with positions n, m of each particle on a chain of length L, hopping probability t (which 
we use to define the energy scale, i.e., t = 1), onsite interaction strength U, and random 
potentials e„ € [—W/2,W/2] for disorder W. Shepelyansky suggested that even for 
repulsive interactions the two particles would form pairs with larger localization length 
A2 than the two separate single particles (SP). Thus the interaction would enhance 
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the possibility of transport through the system. The perhaps even more surprising 
prediction was that at pair energy E = 

\2 cx U^\i^, (2) 

where Ai is the SP localization length. Since Ai oc in ID, this implies large 

values of A2 for small disorders W . 

The first numerical studies devoted to the TIP problem used the transfer ma- 
trix method (TMM) to investigate the proposed enhancement of \2}^'^'^^^^ Other 
direct numerical approaches to the TIP problem have been based on the time evo- 
lution of wave packets, ^'^'^^"^^ exact diagonalization,^^ variants of energy-level statis- 
tics"^*^""^^ and multifractal analysis, '^■''■'^ Green function approaches, pcrturbative 
methods''^'^" and mappings to effective models. ^^^^'"^ In these investigations usually 
(but not always^^"^^) an enhancement of A2 compared to Ai has been found but the 
quantitative results tend to differ both from the analytical prediction in Eq. (2), and 
from each other. Furthermore, a check of the functional dependence of A2 on Ai is 
numerically very expensive since it requires very large system sizes L ^ A2 ^ Ai. Ex- 
tensions of the original arguments have been proposed for higher dimensions, ^^'^^'^^ 
TIP close to a Fermi sea"'^ and long-range interactions in ID^*'^^'^^ and 20.^^^ 

In this paper, we present numerical results for the interaction matrix elements 
and show that a correct statistical description necessitates the use; of the logarithmic 
instead of the commonly used arithmetic average. The resulting typical matrix element 
Utyp has a different disorder dependence than the arithmetic mean Uabs- Following the 
arguments of Ref.,^"^'^'* the use of utyp does not give rise to a power-law enhancement 
as in Eq. (2). However, taking into account the energy denominators, the power-law 
can be recovered but with a smaller power of, e.g., 1.4 ± 0.2 for U = 1. This value is 
in agreement with previous results in m 17-21-24. 27, 28. 30, 31, 3.3-41, 45, 47 -^^c further review 
results based on the application of the decimation mcthod^^ ^'^ which allows for a direct 
computation of \2{U). Augmenting the analysis with the finite-size-scaling (ESS) 
approach, we construct estimates (,2{U) of the TIP localization length in the infinite 
system. ^^'^^ For U ~ they reproduce accurately the well-known dependence of Ai 
on disorder while for finite U, we find ^2{U) ~ 6(0)^^^) with (3{U) varying between 
/3(0) = 1 and /3(1) w 1.5. Thus the enhancement persists, unlike for TMM, 1^^21.23, 24 
in the limit of large system size. Finally, we apply the decimation method in 2D. 
Again using ESS we find scaling functions for 0.5 < U < 2 which exhibit two branches 
corresponding to localized and extended behavior for TIP in an infinite system, and 
we estimate the U dependence of the critical disorder and the critical exponent at the 
transition. 

Before presenting our data, we briefly review the analytical approach to TIP. The 
prediction of Shepelyansky is based upon looking at the interaction matrix element 
of a pair state ipki — i^ki'i with another pair state t/Jnm-^^''^^ Here ipk denotes the SP 
eigenstate localized with Ai around site k and we can restrict the states by — /| < Ai, 
|n — m| < Ai and j/c — n| < Ai, |/ — to| < Ai. The interaction matrix elements are 

JV 

U = ii^kimnm) = Uj24ij)^lU)Mj)i'm{j), (3) 
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Fig. 1 Normalized distributions Po{u) of the off-diagonal coupling matrix elements u for 
Al = 26 {W = 2), L = 200, and U — 1. Thin solid, short-dashed, and long-dashed lines 
correspond to the TIP problem, the TIP problem with random interaction, and Eq. (4), 
respectively. The thick solid line denote TIP matrix elements for states with energy separation 
less than AE = 1/2. Circles, diamonds, triangles, and squares mark the data for the 10 
smallest [u| in each case. The vertical lines above the inset indicate the respective values of 
Wabs > Wtyp as Computed from the total distribution P{u). Inset: same data on a linear scale. 



where we used the Hubbard onsite interaction^ f^Z^j=i ^ji%T with the number op- 
erator rija at site j and spin a. Assuming^^'^^ that the SP state is given as 

Mj) « exp - (^^^ + te,^ (4) 

with dj a random phase, one finds^^ that the average of u has a magnitude of 

Mabs OC A7^^^. (5) 

Shepelyansky calculated the decay rate F of a non-interacting eigenstate by means of 
Fermi's golden rule T ~ U^/Xit}^'^'^'^^ Since the typical hopping distance is of the 
order of Ai the diffusion constant is D ^ U'^Xi/t. Within a time r the particle pair 
visits N ~ C/A^^^f~^/^r^/^ states. Diffusion stops when the level spacing of the visited 
states is of the order of the frequency resolution l/r. This determines the cut-off time 
T* and Aa ~ VDT^ ~ {U/t)'^Xi in agreement with Eq. (2). Applicability of Fermi's 
golden rule requires F ^ t/Xf which is equivalent to U'^Xx/t^ » 1. This is exactly the 
condition for an enhancement of A2 compared to Ai . Alternatively, similar results can 
be obtained by mapping the model onto a random-matrix model with entries chosen 
according to Eq. (5). 13.43,44 
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Fig. 2 Dependence of Wabs (1^) and Mtyp (o) on Ai for the TIP eigenstates from the entire 

band (small symbols) and for AE ~ 1/2 (large symbols) for [7 = 1 and L = 200. The bold 
symbols indicate Mabs and tttyp for L = 250. The solid (broken) lines represent the power 
laws Mabs ~ Ar'-^ (Ar'"^) and wtyp ~ Af^'^^ (Ar''*)- 



2 Numerical results for the interaction matrix elements 

The qualitative arguments presented above should be checked for quantitative accu- 
racy. Even before doing so for Eq. (2), it is already worthwhile to test the validity of 
(5) and the subsequent arguments. It was shown^'' that the assumption of a Gaus- 
sian distribution of u — necessary for taking the arithmetic or r.m.s. average — was 
oversimplified. We calculated^^ the dependence of the averages on Ai and system size. 
To this end, we diagonalized the ID Anderson model for given L and W and com- 
puted u for all suitable states and many disorder configurations. In Fig. 1 we show 
the distribution of u and the values of Utyp and ?iabs- Additionally, wc include results 
for an interaction with randomly varying onsite term U{j) and results for SP states 
with random phase as in Eq. (4). Based on these results, it was shown^^ that due 
to the non-Gaussian distribution of u, one should use the logarithmic rather than the 
arithmetic average for u. However, whereas Uabs oc Af the logarithmic average 
obeys utyp oc Aj~^'^^ as shown in Fig. 2. Following the arguments above, this would 
imply A2 cx X\'^ . We emphasize that this does not mean that there is no enhancement 
of A2. Rather, the results^^ indicate that the arguments of first-order perturbation 
theory^^ capture the physics in a somewhat oversimplified form. Noting the possibly 
quite different energies of the pair states, it appears more suitable to average u only 
over such states which differ in energy by AE < ±\U\/2. This leads to a different Ai 
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Fig. 3 TIP localization lengths ^2 after FSS for C/ = (o), 0.1 (A), 0.2 (O), 0.5 (*) and 1 
(□). The data are for W € [1,6]. The dashed lines show fits according to ^2(11) ~ ^2(0)'', 
the solid line sets the reference for U = 0. Inset: Exponent (3 obtained by a power-law fit to 
the data for U = 0.1, . . . , 1. The filled symbols correspond to the fit (8) with finite c. 



dependence with reduced exponent, e.g., Utyp ^ ^1 ^'^ for C/ = 1, cp. Fig. 2, implying 

A2 0cA^4±°■^ (6) 

in good agreement with previous data based on a wide variety of methods^^'^^"^^'^^'^^'^^'^-'^'^^"*"'^'^^'*^ 
and also with the results of the next section. 



3 Enhancement of the TIP localization length in ID 

The first numerical investigations^^'^^ of A2 used the TMM. Two of us studied the TIP 
problem by a different TMM-'^® at large L and found that (i) the enhancement A2 / Ai 
decreases with increasing L, (ii) the behavior of A2 for [/ = is equal to Ai in the limit 
i — > 00 only, and (iii) for [/ ^ the enhancement A2/A1 also vanishes completely in 
this limit. Therefore we concluded^^"^" that the TMM applied to the TIP problem in 
ID measures an enhancement which is due to the finiteness of the systems considered. 
It is now well understood^'' that the main problem with the TMM approach is that the 
enhanced A2 is expected to appear along the diagonal sites of the TIP Hamiltonian, 
whereas the unsymmetrized TMMs^'^'^* proceed along an SP coordinate. Various new 
TMM techniques have been developed to take this into account, ^^'^'''^^'^^'^^ but still 
all TMMs share a common problem: in general the 11 = result for A2 does not equal 
the value of Ai/2 which is expected for non-interacting particles.^'' Rather, they yield 
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Fig. 4 Schematic view of two possible 
TIP decays in the 2D bar geometry. 



L 



A2(0) w A2(l) and thus A2(0) ^ Ai/2. Therefore other methods appear to be more 
appropriate for the TIP problem. 

Here we briefly review the results''^"^^ of another well-tested method of comput- 
ing localization lengths for disordered systems, the decimation method^^ where one 
replaces the full Hamiltonian (1) by an effective Hamiltonian for the doubly-occupied 
sites only. In accordance with the SP case,^*' A2 is defined by the TIP Green function, 
G2(£;), as35 

^ = -^ln|(l|G.|L)|. (7) 

We computed'*^"''^ the Green function at E = for 26 disorder values W E [0.5,9], 
for 24 system sizes L G [51, 251], and 11 interaction strengths U = 0, 0.1, . . . , 1.0. For 
each triplet of parameters {W,L,U) we averaged I/A2 over 1000 samples. By FSS of 
the X2{U) data we obtained the infinite-size localization lengths £,2{U) shown in Fig. 3. 
We fitted our results by various suggested models. The best fit was obtained with 



1-^ 



6(0) 



(8) 



Such a f7-dependent exponent for the enhancement 6(C^)/C2(0) had been previously 
predicted**^ with /3 up to 2 for U = 1. However, we find that f3 < 1.5 as shown in 
the inset of Fig. 3. Thus we do not sec the quadratic behavior (2) when using the 
fit function of Ref.."*^ On the other hand, after scaling the data onto a single scaling 
curve and fitting 

6(c/)-6(0)cx6(0)'' (9) 

as proposed with (3 — 2,^^ we find^^ indeed /3 = 2 for not too small disorder, e.g., 
W > 2.5 for U = 1, but observe a crossover to /3 = 3/2 for smaller W. 

For U > 1.5 the enhancement decreases again'^'^^^^; the position of the maximum 
depends upon W. We thus find that the duality proposed^'^ for results at U and 
s/^AjU — which would imply a fixed maximum at [/ = 24^^^ — is approximately 
valid. 

We remark that similar results are obtained when placing TIP in different poten- 
tials'*^ which is of relevance for a proposed experimental test of the TIP effect. 
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Fig. 5 FSS scaling curves (lines) and localization lengths A2 for TIP in 2D as a function of 
the scaling parameter ^2 for f/ = (o), 1 (n), and 2 (V). 



4 Localization-delocalization transition for TIP in 2D 



Imry^^'^^ argued that two (onsite-) interacting particles in a 2D random potential show 
enhancement of A2 at = as 



A2 oc Al exp 



Al oc exp 



(10) 



with Al the SP localization length in 2D.^° Indeed, Ortuno and Cuevas^^ numerically 
find a distinct enhancement in 2D for two particles with a short-ranged interaction. 
It is even so strong that an MIT occurs ioY U = 1 at Wc = 9.3 ± 0.2 with critical 
exponent v = 2.4 ± 0.5. Their result is based on the recursive Green function method 
previously employed for the ID TIP case.^^ 

We applied the decimation method to the 2D TIP problem for rectangular bars of 
size M X L with M < L at = for 36 disorders e [3.5, 12], for system widths 
M = 2, . . . , 8 at fixed L = 51, and 51 interactions strengths U = 0, 0.04, . . . , 2.0 with 
hard wall boundary conditions. For the results presented here, we average the decay 
length of TIP states over every possible transmission path between coordinates (1,1) 
and various heights y on the right hand side (L, y) of the bar as shown in Fig. 4 



1 

A^ 



1 ^ 
u 51 



M 



lnKl,l|G2|X,j/)| 

^ V(i-l)2 + (j;-l)2 



(11) 



w 



Here {■)w denotes the disorder average over 1500,500,250,200,100,50 and 20 con- 
figurations for M = 2,..., 8, respectively.^^ To obtain a more accurate result the 
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Fig. 6 Phase diagram of the TIP localization-delocalization transition in 2D. The data point 
(□) for ?7 = 1 indicates the result of Ref..^* The error bars represent the error propagation 
of the Levenberg-Marquardt nonUnear fit. Inset: power-law fit to the data for U > 0.4 with 
exponent 0.36. 



value of I/A2 was actually obtained by fitting the slope of (In |(1, l|G2|a;,t/)|) versus 
a/(x - 1)2 + {y- 1)2 where x < L. 

In Fig. 5, we show typical FSS results for three values of U. The quality of the 
scaling curves is not as good as in the ID TIP analysis, due to the smaller samples 
and smaller number of configurations. The scaling curves for U < 0.2 can be described 
by a single branch corresponding to localized behavior. For U > 0.5, the scaling curves 
have two branches indicating a transition of TIP states from localized to delocalizcd 
behavior. In the intermediate range 0.2 < U < 0.5, FSS is less convincing and requires 
further investigations. We remark that the FSS results do not change appreciably 
within the accuracy of the data when we drop the data for A'/ = 2 or 8. 

In order to reliably extract the critical disorders Wc and the correlation length 
exponents from the FSS data, we employed the FSS techniques of Refs.^° and.^^ In 
the lowest order approximation, the fit function is given as^° 

^=A.+.«./«(|;-i), (12) 

with fit parameters Ac, a, Wc and f. We used a fit function containing up to the 
third power in W/Wc-^^ In this way including deviations from the linear behavior in 
W is important due to limited quality of the data. In Figs. 6 and 7, we show the 
resulting dependence of Wc and v on U. Note how the FSS gives unreliable results 
for U < 0.5, resulting in the incorrect prediction of Wc(0) > 0. As demonstrated in 
the inset of Fig. 6 the Wc data can be approximated for 0.4 < ?7 < 1 by a power-law 
fit Wc oc [/0-36±o.o3 rpjjg critical exponent p shows a slight decrease for U > 0.5 and 
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Fig. 7 Critical exponent v according to Eq. (12). The data point (□) for f/ = 1 represents 
the result of Ref..^* The error bars are obtained as in Fig. 6. 

is in reasonable agreement with the value of 2.4 ± 0.5 obtaincxi in Ref.^'' for [7 = 1. 
However, the (non-universal) value of Wc = 8.1 ± 0.1 at i7 =1 is somewhat different 
Ref..^^ We attribute this to our different definition (11) of A2.^^ 

5 Conclusions 

In conclusion, we observe an enhancement of the TIP localization lengths due to onsite 
interaction both in ID and 2D random potentials. This enhancement persists, unlike 
for TMM,^^'^^'^"^ in the limit of large system size and after constructing infinite-sample- 
size estimates from the FSS curves. In the 2D case, it leads to numerical evidence for 
a localization-delocalization transition of the TIP states. We emphasize that this is 
not a metal-insulator transition in the standard sense since only the TIP states show 
the dclocalization transition. The majority of non-paired states remains localized. 

RAR gratefully acknowledges support by the DFG via SFB-393. ML gratefully acknowledges 
the support of the EU-TMR programme under grant no. ERBFMBICT972832. 
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